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We use a magnetic force microscope (MFM) to investigate single vortex pinning and penetration 
depth in NdFeAsOi-ccFcc, one of the highest-Tc iron-based superconductors. In fields up to 20 Gauss, 
we observe a disordered vortex arrangement, implying that the pinning forces are stronger than the 
vortex-vortex interactions. We measure the typical force to depin a single vortex, Fdepin — 4.5 pN, 
corresponding to a critical current up to Jc — 7 x 10^ A/cwS. Furthermore, our MFM measurements 
allow the first local and absolute determination of the superconducting in-plane penetration depth 
in NdFeAsOi-ccFa;, \ah = 320 zb 60 nm, which is larger than previous bulk measurements. 


PACS numbers: 74.25.Wx, 74.70.Xa, 68.37.Rt 

I. INTRODUCTION 

A long-standing challenge on the road to supercon¬ 
ducting applications is the unwanted motion of vortices, 
quanta of magnetic flux 4>o = h/2e which penetrate a 
type-II superconductor in the presence of a magnetic 
field. Each vortex consists of a circulating supercurrent 
which decays radially on the length scale of the penetra¬ 
tion depth A, and a core of suppressed superconductiv¬ 
ity with radius given by the coherence length ^o- When 
a current is applied to the superconductor, the vortices 
experience a Lorentz force density Fp = J x B, perpen¬ 
dicular to the direction of the current and proportional 
to its magnitude. Although the net supercurrent remains 
non-dissipative, the accompanying motion of the normal 
electrons in the vortex cores causes resistance and energy 
loss in the material. 

A vortex can be pinned by co-locating its energetically 
costly core with a preexisting crystal defect where super¬ 
conductivity is already suppressed. The strength of the 
pinning force density Fp dictates the maximum super¬ 
current which can be applied without vortex motion and 
consequent dissipation. Several decades of engineering 
effort have been devoted to optimizing vortex pinning 
in superconductors.However, vortex pinning in the 
highest-Tc cuprate superconductors remains challenging, 
due in part to the large electronic anisotropy which al¬ 
lows vortices to bend and depin in pancake fragments 
rather than as one-dimensional semi-rigid objects.^ 

The recent discovery of high-Tc iron-based supercon¬ 
ductors (Ee-SCs) brought new optimism to the vortex 
pinning problem.Ee-SCs were found to be typically 
more isotropic than their cuprate cousins.Eurther- 
more, intrinsic pinning tests showed promise,^ raising 
hopes that defect engineering could improve the pin¬ 
ning properties to achieve critical currents larger than 
the cuprate benchmark around 10^ A/cm^.^ 

Maximizing critical current will require detailed under¬ 


standing of the pinning efficacy of specific defects. Trans¬ 
port measurements yield information about macroscopic 
critical currents but cannot directly distinguish the dis¬ 
tribution of pinning forces and the sites responsible for 
the strongest pinning. Eurthermore, the dissipative mo¬ 
tion of small numbers of vortices may not be detected by 
macroscopic transport measurements. It is crucial to un¬ 
derstand and quantify pinning at the single vortex level. 
An important parameter in this endeavor is A, the funda¬ 
mental length scale of magnetic interactions, whose abso¬ 
lute value has remained challenging to measure via bulk 
techniques due to pervasive inhomogeneity in Ee-SCs. 

While scanning tunneling microscopy, Bitter 

decoration,and scanning techniques 

can image vortex configurations and offer some insight 
into pinning, magnetic force microscopy (MEM) is the 
most powerful technique to directly quantify the single¬ 
vortex depinning force Tdepin and the in-plane penetra¬ 
tion depth Xab- MEM was used to measure vortex lat¬ 
tice correlation length, Tdepin, and Xab distributions in 
doped BaEe 2 As 2 ,^^~^^ but these parameters have not 
been locally quantified in the highest-Tc ‘1111’ family 
of Ee-SCs. More details of previous work are summa¬ 
rized in Appendix A. Here we use MEM to investigate 
NdEeAsOi-a^Ea^. We find a typical single-vortex de¬ 
pinning force Tdepin — 4.5 pN and penetration depth 
Xab = 320 ± 60 nm. The tip-induced vortex motion 
demonstrates out-of-plane electronic anisotropy, but no 
statistically significant evidence of in-plane anisotropy in 
NdEeAsOi-a^Ea,. 


II. EXPERIMENT 

The single crystal used in this experiment, of thick¬ 
ness ^10 fim and lateral dimension ^100 /im, shown 
in Eig. 1(a), was mechanically extracted from polycrys- 
talline NdEeAsOi-a^E^, synthesized at high pressure. 
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FIG. 1. Sample characterization, (a) Photograph of the 
NdFeAsOi-ccFa; single crystal, (b) Tunnel diode resonator 
experiment shows Tc(onset) ~ 50K. (c) SEM image of the 
same crystal, (d) Room temperature EBSD Kikuchi patterns 
determining the axis orientation. This sample is expected 
to remain in the tetragonal structure for all temperatures,^^ 
but our EBSD data lacks the resolution to independently 
distinguish between tetragonal (P4/n?nm) and orthorhombic 
(Cmma) structures. 


It is slightly underdoped with nominal x = 0.1, and 
Tc(onset) 50K measured by tunnel diode resonator 
(Fig. 1(b)). The sample was cleaved in air perpendicular 
to the c-axis, aligned in the MFM, pumped to high vac¬ 
uum, cooled to base temperature ^6 K, and imaged in 
cryogenic ultra-high vacuum (UHV) at applied fields up 
to ±20 Gauss. After the MFM experiment, the sample 
crystallinity and orientation were characterized by scan¬ 
ning electron microscopy (SEM) and electron backscatter 
diffraction (EBSD) as shown in Figs, l(c-d). 

MFM imaging was carried out in frequency modula¬ 
tion mode^^ using a commercial NSC18 cantilever from 
Mikromasch with manufacturer-specified force constant 
k = 3.5±1.5N/m.^^ As the tip oscillates in the ^ axis, the 
magnetic moment at the end of the tip interacts with the 
vortex. The bare cantilever resonance frequency /o = 73 
kHz is shifted by the local force gradient according to 

fodF^{x,y,z) 

^ 2k dz ' ^ ^ 

The total force consists primarily of van der Waals and 
magnetic components, with the former dominant when 
the tip is close to the sample.We use the rapid change 
in dFz / dz to determine the position of the sample surface 
to within a few nanometers, then retract the tip by a fixed 
amount 2 ; for constant height imaging in a regime where 
the magnetic force is dominant. 

We employ a feedback loop to maintain the tip oscil¬ 
lation at its local resonance frequency as the sample is 



FIG. 2. MFM images of vortices vs. applied field Ba in 
NdFeAsOi-ccFcc. (a) Ba = —9 G {z = 200 nm). (b) Ba = —4 
G (z = 300 nm). (c) Ba = I G {z = 400 nm). (d) Ha = 3 G 
{z = 300 nm). (e) = 8 G (z = 300 nm). (f) Vortex density 

vs. Ba- At these scan heights, the tip-vortex force does not 
cause significant vortex depinning. 


scanned; the resultant map of the local frequency shift 
df thus serves as a map of the vertical force gradient, 
dFzjdz. Fz is not directly imaged, but it can be obtained 
by integrating the measured dFzjdz. The lateral force 
components can then be estimated by assuming a trun¬ 
cated cone tip shape.Although Newton’s third law, 
which dictates that any force used as an imaging sig¬ 
nal also influences the sample itself, is typically regarded 
as a disadvantage which makes a force microscope an 
invasive probe, we use it here to our advantage in or¬ 
der to pull or push on vortices and directly measure the 
forces required to dislodge them from their pinning sites 
in NdFeAsOi_a,Fa,. More details of the tip-vortex inter¬ 
action regimes are given in Appendix B. 


III. RESULTS 
A. Vortex Imaging 

Figs. 2(a-e) show images of individual vortices vs. ap¬ 
plied field Ba in NdFeAsOi-a^F^,. Before the acquisition 
of each image, the sample was heated above Tc, to 60 
K, then cooled to 6 K within the field shown. We note 
immediately that the vortices are disordered, indicating 
that at this low density the pinning forces on individ¬ 
ual vortices dominate over the vortex-vortex interaction 
forces in this NdFeAsOi_a,Fa, sample. The vortex den¬ 
sity vs. applied field is shown in Fig. 2(f). The intercept 
serves to calibrate the ambient B field offset, while the 
slope serves to calibrate the (x, y) scan piezo. We And an 
offset field of 1.8 Gauss, consistent with Earth’s field and 
stray fields from steel screws close to the MFM. For all 
following figures and discussion, we employ the calibrated 
B and (x, y) distances. 









3 


(a) z = 350nm (b) z = 300nm (c) z = 250nm (d) z = 200nm (e) z = 150nm 
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FIG. 3. Vortex depinning. (a-e) MFM images of the same field of view containing ~9 vortices at T = 6 K and B = 5.8 Gauss, 
with decreasing tip-sample separation z. As decreases, the increasing lateral force between tip and sample causes some 
vortices to depin. Black circles indicate vortices depinning for the first time, (f) Maximum vertical force gradient {dFz/dz)max 
vs. z + Xab (with Xab = 320 nm) for the 7 vortices entirely within the held of view. The values for each vortex are depicted 
as a different series of colored dots, terminating at the minimum z where the vortex remains pinned. Blue line shows a ht to 
Eq. 2 with hxed ho = 200nm and free parameters A and Xab- The gray shading denotes the region between {z + Xab)~^ and 
{z + Xab)~‘^: which is used to bound any systematic error in the truncated cone tip model. Green dashed lines indicate the 
other major source of systematic error, the uncertainty in the manufacturer-specihed spring constant k. Inset shows a sketch 
of the magnetic tip, which dehnes ho. (g) The integral of (f), using the functional form of the ht to extrapolate to z = oo, gives 
the total shown as a blue line. Gray shading and green dashed lines dehne the systematic errors from the model and the 

spring constant, respectively, (h) Histogram of Fhepin for the 6 displaced vortices shown in (a-e). 


B. Pinning Force Distribution 


We measure the vortex pinning force distribution by 
acquiring a series of images starting at tip height = 
950 nm, and approaching the sample, as exemplihed in 
Figs.3(a-e). As the tip-vortex force increases in succes¬ 
sive scans, we note the height ^^depin of the hrst depinning 
event for each vortex (marked with a circle). To obtain 
the force corresponding to ^^depin, we hrst plot in Fig. 3(f) 
the maximum force gradient {dFz/dz)^^,^ at each vortex 
center. To compute Fz at these locations, we must inte¬ 
grate dFzjdz from ^ = oo to ^^depin- Between z^^x and 
Zpin, we simply sum the measured {dFz/dz)max- To han¬ 
dle the tail from 2 ; = oc to z^a^^ we model the vortex as a 
monopole situated Xab beneath the surface.We model 
the tip as a uniformly coated truncated cone, depicted as 
an inset to Fig. 3(f),^^ which leads to 


dz 


= ^oA 


2ho 


_ + Xab^ + Xab)^ _ 


( 2 ) 


Here ho is the tip truncation height (which we hx at 200 
nm based on manufacturer specihcations) and A = atm 
is a tip-specihc ht parameter depending on the opening 
angle a, the magnetic him thickness t, and the magneti¬ 
zation m. We ht the measured {dFz/dz)max vs. 2 ;, extrap¬ 
olate the ht function to 2 : = 00 , integrate the tail, and 
add this result to the sum of the measured (dFz/dz)max, 


to arrive at the total shown in Fig. 3(g). We take 

Tdepin to be the maximum lateral force which is 

a factor of ^ 0.31 times F^^^ for a truncated cone tip. 
Finally, Fig. 3(h) shows a histogram of Fdepin for the six 
vortex motion events within this held of view. Details of 
the analysis and error bars are given in Appendix C. 

From Fig. 3(h), the typical single-vortex depinning 
force is Fdepin — 4.5 pN. In order to convert the pinning 
force to a critical current Jc = Fdepin/(4>o^), we need an 
estimate of the length £ of the vortex which is depinned. 
Although the sample thickness is on order 10 /im, which 
places an upper bound on our observation of vortex 
wiggling in Fig. 5(b), as well as the occasional observa¬ 
tion of a hux ‘tail’ after a vortex dragging event in Fig. 
5(i), suggests the possibility that only the top portion of 
the vortex is depinned by the tip. The visibility of the 
stationary bottom portion of the vortex decays exponen¬ 
tially with its distance beneath the surface,placing 
an approximate lower bound of A on the length of the 
upper displaced portion. The measured Fdepin — 4.5 pN 
thus corresponds to 4.5 x 10“^ to 1.4 x 10“^ N/m. 

The measured depinning force is roughly consistent 
with the density ^2-60x10“^^ m“^ of strong pinners 
with individual strength 2 x 10“^^ N per pinner previ¬ 
ously estimated by magneto-optical imaging on under¬ 
doped NdFeAsOi-a^Fa, (Tc ^ 35 K).^^ Our measurement 
therefore implies a critical current 2 x 10^ A/cm^ ^ Jc ^ 
7 X 10^ A/cm^. These values greatly exceed the bulk crit- 

















4 


ical current density Jc ^ 2000 A/cm^ in polycrystalline 
NdFeAsOi-a^Fa^,^^ and are roughly consistent with the 
in-grain persistent current density detected by magneto¬ 
optical imaging, Jc r\j 10^ A/cm^ (Refs. 43 and 41) and 
Jc ^ 5 X 10^ A/cm^ (Ref. 42), and with requirements for 
various superconducting applications.^ 


C. Penetration Depth 


The superconducting penetration depth A is necessary 
to quantitatively model the single vortex pinning force, 
as well as to determine fundamental parameters such as 
the phase stiffness and superfluid density. The tempera¬ 
ture dependence of A(T) can be used to distinguish be¬ 
tween fully-gapped vs. nodal pairing scenarios.How¬ 
ever, the absolute value of A has been notoriously diffi¬ 
cult to measure in 1111 Fe-SCs via bulk techniques, due 
to prevalent sample inhomogeneity. 

Two common techniques to measure A are muon spin 
rotation (/iSR) and lower critical field (i^ci) measure¬ 
ments. First, /iSR experiments measure a histogram of 
the spin decay times of stopped muons within a supercon¬ 
ducting sample, which translates to a histogram of the 
local magnetic field strengths present in the sample. As¬ 
suming a triangular vortex lattice in NdFeAsOo.ss with 
Tc = 51 K, one /iSR study extracted Xab = 195 ±5 nm,^^ 
but our observation of a disordered vortex arrangement 
in NdFeAsOi-a^Fa, suggests that this /iSR study under¬ 
estimates Xab- Second, Xab can be related to via 
Hci = ^o/(47rA^5) [In/i: + c(/i:)], where n = Xab/^ab and 
c{hz) is a function that tends to 0.5 for large values of A 
measurement of = 150 Gauss and ^ 2.9 nm (from 
Hc 2 = 40 T) can be used to determine Xab = 270 ± 40 
nm in NdFeAsOi-a^F^, with Tc = 32.5 

We use MFM to measure Xab locally, via the Meissner 
repulsion of the tip in a vortex-free region. For z > Aq^, 
the tip-superconductor interaction can be approximated 
by the interaction between the tip and an image tip re¬ 
flected across a mirror plane a distance Xab below the 
surface of the superconductor.^^ For z > 300 nm, we fit 
to the Meissner repulsion 


OF, OF, 
dz dz 

= 2'kijlqA^ 


z=oo 

1 

^ + ^ab 


+ 


ho hi 

{z + XabY 2(2: -h Xab)^ 


( 3 ) 


to find Xab = 320 ± 60 nm. Details of the analysis are 
given in Appendix D. 

From our absolute measurement of A^^, we can com¬ 
pute several fundamental quantities. Using ksTp = 
h?'KC\atnslrrC and 1/A^^ = jc? • where 

Qnt = 8.557 A is the distance between superconduct¬ 
ing layers,we compute = 650 K, which is in rough 
agreement with the trend Tp ^ lOTc noted for a variety 
of unconventional superconductors.^^ 



FIG. 4. Penetration depth, (a) The blue line shows df vs. 
tip-sample distance z at a location which is laterally at least 
2 /xm from the nearest vortex (T = 6 K; R = 1 Gauss.) The 
red line shows a fit to Eq. 3 for z > 300 nm. (b) Blue dots 
show the background values of df from the vortex images in 
Fig. 3. The red line is a fit to Eq. 3 for 2 ; > 300 nm. 


IV. CONCLUSION 

In summary, we have presented a systematic inves¬ 
tigation of single vortex pinning in one of the highest- 
Tc Fe-SCs, NdFeAsOi-a^Fa,, for B up to 20 Gauss. We 
find a disordered vortex arrangement, implying that the 
pinning forces are stronger than the vortex-vortex in¬ 
teractions. The average pinning force of a single vor¬ 
tex is Fp 4.5 pN, consistent with a critical cur¬ 
rent jc up to 7 X 10^ A/cm^, within striking distance 
of the benchmark value for cuprates. NdFeAsOi-a^F^, 
presents a rich opportunity for vortex physics, with the 
out-of-plane anisotropy sitting right at the boundary be¬ 
tween pancake and rigid Abrikosov vortices. (No in-plane 
anisotropy was detected, as detailed in Appendix E.) 
NdFeAsOi-a^Fa, also presents significant technical oppor¬ 
tunity to improve the inter-grain jc by the reduction of 
wetting impurity phases^^ and the intra-grain jc by the 
improvement of nanoscale pinning defects. The ability to 
intentionally position the vortex top using the MFM tip 
opens the door for experiments to measure pinning forces 
at specific sample defects which may be independently 
imageable in AFM mode (i.e. with smaller tip-sample 
separation 2 ;, where topographic and electrostatic forces 
dominate over magnetic ones). 
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Appendix A: Previous Vortex Imaging Studies 

Single vortices have been directly imaged in all four 
major families of Ee-SCs. Due to the relative ease of large 
single crystal growth, the ‘ 122 ’ family has been the most 
heavily investigated. Scanning tunneling microscopy 
(STM) first showed a disordered vortex lattice indicative 
of strong pinning in optimally doped Ba(Eeo. 9 Coo.i) 2 As 2 
with Tc = 25.3 Eurthermore, the lack of correla¬ 
tion between vortex locations and surface defects sug¬ 
gested that vortices behave as semi-rigid one-dimensional 
objects,consistent with low anisotropy 7 = ^ 

2 measured by transport techniques.^ Bitter decoration 
imaging of larger arrays of vortices in the same ma¬ 
terial confirmed the vortex disorder,while mag¬ 
netic force microscopy (MEM) quantified the scaling 
of vortex disorder with applied field.In underdoped 
Ba(Eeo. 95 Coo.o 5 ) 2 As 2 (Tc = 18.5 K), MEM imaging 
demonstrated the uniformity of the superfluid density at 
the length scale of the penetration depth \ab = 325 ± 50 
nm, and allowed measurement of a typical single-vortex 
depinning force Tdepin = 18 pN at T = 5 K in a twin- 
free, 10 /im-thick sample, corresponding to 1.8 x 10“^ 
N/m.^^ Another Bitter decoration study of under- and 
optimally-doped Ba(Eei_a^Coa^) 2 As 2 found an even larger 
typical pinning force of order 10“^ N/m.^^ In contrast, 
STM imaging of optimally doped Ko. 4 Bao. 6 Fe 2 As 2 with 
Tc = 38 K showed an ordered vortex lattice indicative 
of weak pinning.MEM imaging of optimal and under¬ 
doped Ka,Bai_a,Ee 2 As 2 showed disorder at low applied 
fields, with a tendency towards an ordered lattice at Tesla 
fields.Eurthermore, the low-field (10 Gauss) pinning 
force was measured to be 2 x 10“^ N/m,^^ a full order 
of magnitude weaker than in Ba(Eeo. 95 Coo.o 5 ) 2 As 2 .^^ A 
possible resolution to the pinning differences within the 
‘122’ family was prompted by the observation of disor¬ 
dered vortices Ka,Sri_a,Ee 2 As 2 :^^’^^ dopant size mismatch 
and consequent clustering may increase electronic inho¬ 
mogeneity at the ^0 length scale and enhance vortex 
pinning.Very disordered vortices were also observed 
by Bitter decoration in Ba(Eei_a^NU) 2 As 2 , with a ten¬ 
dency to form vortex stripes in some regions.Einally, 
MEM imaging of vortices in Ba(Asi_a;Pa^) 2 Ee 2 showed 
relatively strong, uniform pinning for overdoped samples, 
clustered pinning along twin boundaries for slightly un¬ 


derdoped samples, and weak pinning for very underdoped 
samples. 

In the ‘11’ family, STM imaging showed a disordered 
vortex arrangement in EeSca^Tei-a^ with x ^ 0.4,^^ but 
an ordered vortex lattice in stoichiometric EeSe films, 
consistent with the idea that dopant clustering may cause 
the inhomogeneity necessary for pinning. 

In the ‘ 111 ’ family, STM imaging of LiEeAs found 
vortex disorder increasing with applied field,which 
can be attributed to a decreasing shear modulus Cqq 
that allows pinning forces to dominate over the vortex- 
vortex interaction forces.STM imaging of vortices in 
Na(Eeo. 975 Coo.o 25 )As at 11 T also showed a slightly dis¬ 
ordered lattice. 

Scanning squid microscopy images of vortices in 
twinned Ba(Eei_a;Coa;) 2 As 2 showed that vortices avoid 
twin boundaries,consistent with a report that twin 
boundaries enhance the local superfluid density.In con¬ 
trast, twin boundaries were found to pin vortices in 
EeSe.^^ Taken together, all of these studies show diverse 
pinning behavior across the ‘ 122 ’, ‘ 11 ’, and ‘ 111 ’ families 
of Ee-SCs. 

In the ‘ 1111 ’ family, scanning SQUID mi¬ 
croscopy showed disordered vortices in polycrystalline 
NdEeAsOi-a^Ea,, but the vortex profiles were severely 
resolution-limited.^^ Bitter decoration of single crys¬ 
talline SmEeAsOi-a^Ea, at 46 Gauss also showed 
disordered vortices, with some tendency to form stripes, 
but the images were again severely resolution-limited.^^ 
Neither the intrinsic vortex profile, nor the single 
vortex pinning forces have ever been quantified in this 
highest-Tc famiiy of Ee-SCs. 


Appendix B: Tip-Vortex Interaction Regimes 

We tune the tip-vortex force by varying the tip height 
z, as shown in Eig. 5. At iarge z (‘surveiiiance height’), 
the interaction force is smaii, so the MEM can be used 
to image vortices without dispiacing them, as shown in 
Eigs. 5(a-c). 

As we iower the tip height to a regime which exceeds 
most pinning forces, we observe that most vortices wiggle 
back and forth as the tip is rastered over them, as shown 
in Eigs. 5(d-f). This behavior implies that the vortex is 
bending at the top, as sketched in Eig. 5(d). 

As we lower the tip still further, turning off the oscilla¬ 
tion and approaching to within 50 nm of the sample, the 
tip-sample force becomes strong enough to interact with 
a deeper portion of the vortex, and the full observable 
top of the vortex can be moved as a single unit, as shown 
in Eig. 5(g-i). However, we do occasionally see weak tails 
which indicate that even in this interaction mode, the 
top portion of the vortex may decouple from the lower 
portion, suggesting significant out-of-plane anisotropy. 

The energy cost of an Abrikosov vortex scales with the 
volume of its core and the free energy of the superconduc¬ 
tor. To minimize the core volume, the vortex will take the 
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FIG. 5. (a-c) Large z: passive observation, (d-f) Intermediate z: only the top of each vortex is displaced, as it stretches away 
from an anchor point deeper within the crystal, (g-i) Small z: in some cases, the entire observable vortex is displaced. Note 
that one of the vortices has a ‘tail’ indicating that only the top portion of length ~ A was displaced with respect to the bottom 
portion. 


shortest path through the sample, deviating only slightly 
to take advantage of pinning sites along the way. How¬ 
ever in an anisotropic layered material a vortex may kink 
its way along an indirect path through the sample like an 
offset stack of pancakes.^ Each pancake can pin indepen¬ 
dently, and the tip interacts with each pancake with ex¬ 
ponentially decaying force.In NdFeAsOi_a,Fa, with 
nominal x = 0.18 and Tc(onset) = 50 K, transport 
measurements on single crystals in field yield anisotropy 
7 = H^ 2 l^c 2 — 4-34 — 4.9 while cleaner crystals with 
Tc(onset) = 52 K yield 7 < 6 (Ref. 7). On another sam¬ 
ple of NdFeAsOi-a^Fa, with Tc ^ 34 K, the anisotropy 
was 7 ^ 5.5 from the ratio of the upper critical fields, 
7^4 from the ratio of the lower critical fields,and 
7 ^ 7.5 from the ratio of the irreversibility line.^^ 


Appendix C: Pinning Force Analysis 


The magnetic field of a vortex at a distance z above 
the ab surface of a superconducting half-space 


B{R, z) 


4)0 f o 


(Cl) 


where Xab is the in-plane penetration depth and Q = 
y/l + {Xab(lY‘ The total interaction energy between an 


MFM tip and a vortex is 

U{R,z)= j dz'<fiR'M^{R',z')B^{R + R',z + z') 

J tip 

(C 2 ) 

where Mz{R'^z') is the magnetic density distribution, 
and the primed coordinates refer to the tip, with origin 
at the center of the plane which truncates the conical tip. 
The signal measured by the MFM is the vertical gradient 
of the magnetic force, which is given by 

BTC r)^ 

-^{R,z)= / dz'(fR'M,{R',z')^B,{R+R' ,z+z'). 

(C3) 

For the truncated cone tip model with small opening 
angle a, Eq. C3 can be approximated by 


dz 


{R,z) 


$ 0 ^ 


2 : + ^ab 

.((^ + A„b)2 + i?2)3/2 


ho{2{z + Xabf - R^)] 

((^ + Aa6)2 + i?2)5/2 

(C4) 


where A = atm. Setting i? = 0 in equation (C4) gives 
the peak vertical force gradient above a vortex, 


dF, 


dz 


max 


^oA 


1 

_{z + Xab)"^ 


2ho 

(z + AQb)3 


(C5) 
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FIG. 6. Extraction of {dFz/dz)ma.^ for a single vortex, (a) 
Vortex image (raw data minus plane background) acquired at 
tip height = 450 nm. (b) Fit of truncated cone model in 
Eq. C4 to data in (a), (c) Contribution of a single vortex, 
obtained by subtracting the fits of all-but-one vortex in (b) 
from the data in (a), (d) Cut through the single vortex in 
(c), along the dashed line. Red dots are data, blue line is 
smoothed data, and (dFk/dz)max is taken to be the maximum 
of the blue curve. 



FIG. 7. Method to determine vortex movement, (a) A vortex 
image acquired at tip height z = 300 nm. (b) The same vortex 
at tip height z = 250 nm. (a) and (b) are plotted on the same 
color scale, (c) The difference between the vortex images in 
(a) and (b). (d) Log-log plot of the maximum pixel value of 
each successive difference image, vs. height z. Blue points are 
extracted from the height series data. Red line indicates the 
theoretical upper bound determined from Eq. C4, with up to 
one pixel lateral shift, as described in the text. Arrow shows 
the point corresponding to z = 250 nm, where the slope of 
the data (blue) exceeds the theoretical slope (red). 


We first fit each image in the height series (exemplified 
in Figs. 3(a-e)) to Eq. C4. We fix = 200 nm according 
to tip manufacturer specifications, and fix ^ according 
to each measurement. The fit parameters are the pene¬ 
tration depth Aa 5 , the tip parameter A = atm^ and the 
center coordinates {xi^yi) of each vortex. Example raw 
data and fit are shown in Figs. 6(a-b). We then subtract 
the fit image of all vortices but one, to leave a single vor¬ 
tex, exemplified in Fig. 6(c). This procedure effectively 
removes the tails of the surrounding vortices so that the 
remaining single vortex can be analyzed individually with 
higher accuracy. From this individual vortex, we read off 
the value of {dFz/dz)ma.^^ as demonstrated in Fig. 6(d). 

The values of {dFz/dz)ma.^ for seven vortices fully 
within the field of view are plotted as a function of 2 : 
in Fig. 3(f). We then fit all of these points to Eq. C5, 
with fixed Hq = 200 nm and fit parameters A and Aab- 
This allows us to extrapolate the force gradient to infin¬ 
ity. We integrate this curve from ^ = 00 to the maximum 
measured tip height Zmax, then sum the measured values 
of (dFz/dz)max, to arrive at the total vertical force ex¬ 
erted by the vortex on the tip above the center of the 
vortex, shown in Fig. 3(g). 

To find the depinning forces, we assume that the vor¬ 
tex moves when the tip is scanned over the radius R 


where the lateral force is greatest. The analytic form of 
the lateral force is computed by taking the lateral deriva¬ 
tive of Eq. C2. The ratio of F^^^/F^^^ obtained from 
the analytic form is ^0.31 for ho and 2 ; in our range of 
interest. 

To determine the height at which each vortex is de- 
pinned, we look at the differences between the images 
taken at successive tip heights. This shows any abrupt 
change in vortex shape which indicates vortex movement. 
To make a quantitative determination, we use the vor¬ 
tex model in Eq. C4 to compute the maximum theoreti¬ 
cal value of each difference image between two successive 
heights, allowing also for the possibility that successive 
images are misaligned by one pixel to account for possi¬ 
ble tip drift between images. We plot the measured data 
(blue dots) and theoretical values (red curve) on a log- 
log plot to highlight their difference. A measured slope 
of absolute value greater than the theoretical slope in¬ 
dicates an abrupt change larger than can be accounted 
for by tip height and image misalignment alone. This 
method is depicted in Fig. 7. Using this method, we 
obtain the histogram in Fig. 3(h). By our force curve 
above, we obtain an average vertical force of 14.5 ± 1.2 
pN at the center of the vortex, and an average maximal 
lateral force of 4.5 ± 0.3 pN. Including the systematic 
uncertainty of the cantilever spring constant, this value 



















FIG. 8. Successive approximations of Eq. D5 for the Meissner 
response in a truncated cone tip model. Parameters plotted 
here are ho = 200 nm, ji^A — 1 /xm/Hz, and \ah — 300 nm. 
(We repeated the convergence test for values of \ah up to 350 
nm; in each case the series was found to converge well for 
^ ^ Aab.) 


becomes 4.5 ± 2.5 pN. Uncertainty in the tip truncation 
height ho contributes insignificantly to the total error. 

We note that the depinning forces we report are an 
upper bound on the actual depinning forces, since we 
measure only at 50 nm 2 ; increments, and we assume that 
each vortex moves where the lateral force from the tip is 
greatest. 


Appendix D: Penetration Depth Analysis 


We model the Meissner response following Ref. 54, 
which we reproduce here for completeness. The Meissner 
response of a type-II superconductor exerts a force on an 


MFM tip above the ab plane and far from any vortices, 
given by^^ 


F{z) = ^ [ dkk^G{\abk)e-'^^^ f dr'x 

Jo Jtip 

f dr"M{r')M{r")e~'^^^'+^'"^Mk\R' - R''\) (Dl) 
J tip 

The MFM tip measures the vertical force gradient 

/ dkk^G{Xabk)e~'^^'' f dr'x 
OZ Z-K Jq Jtip 

f dr"M{r')M{r")e-’^^^'+'^"'>Jo{k\R' - R"\) (D2) 

J tip 

where G{x) = (Vl + x‘^ — x)/{x/l + x‘^ + x). Using the 
approximation Jo(h|i^ — i^'|) ^ Jo(0) and 


x/l-\-x‘^ - X 

G(x) = , -^ e 

Vl + + X 


/y*3 

- 2 ^ I l + ^_^ + Olxf 

3 20 ^ ^ 


for X ^ 1, we obtain, to first order G{x) ^ e 


-2a: 


(D3) 


= -27r/io^^x 


1 


+ 


ho 


hi 


Z + Xab + ^ab)‘^ + '^ab)^ 


(D4) 


where A = atm is the tip parameter. 

To confirm the convergence of the series, we also com¬ 
pute the expression including higher order terms in the 
expansion of G{x) from Eq. D3. 


^^( 2 ;) = -2ttiioA^ 


1 


+ 


ho 


+ 


hi 


+ 


{z + \ab) {z + XabY 

1 ho 


Mb 


{z + XabY 

{Z + Aab)^ 


+ 




bhl 


4 (^ + \ab) (z + ^ab)^ ^(z + ^ab)^ 

9 27ho , 189hl 


+ 


+ 


16(^ + Xab) 8{z + Aab)^ 32 ( 2 : + Xab)^ 


+ 0 




{z + Xab)^ 


(D5) 


We plot several partial sums for test Xab and ho values, 
exemplified in Fig. 8, which shows that the series con¬ 
verges rapidly in the range 2 ; > A^b, while it does not 
converge well for 2 ; ^ Xab- Thus we are justified in using 
the first order approximation (Eq. D4) for 2 ; > Xab- 

To measure the Meissner response curve, we scan the 
tip from high to low above the sample at a fixed location 
^2 /xm laterally from the nearest vortex, as plotted in 
Fig. 4(a). We see that the force curve changes abruptly 
for 2 ; < 100 nm. We attribute this feature to electrostatic 
forces between the tip and the sample. However, this 


force gradient falls off as 1 , so we may neglect it for our 

tip heights of interest 2 ; > 200 nm, where the dominant 
Meissner force gradients fall off as Ijz. 

We fit the Meissner response curve in Fig. 4(a) from a 
cutoff height of Zcut = 300 nm to Eq. D4, with fixed ho, 
and fit parameters A and Xab- We find Xab = 320T60 nm. 
The greatest uncertainty comes from our uncertainty in 
the tip truncation height ho, which is nominally 200 nm 
according to manufacturer specifications, but which we 
have not measured directly. Thus the error bars are 
determined by redoing the fit with ho = 100 nm and 
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FIG. 9. Vortex wiggling anisotropy, (a-d) Series of four rotated images exemplifying the data acquisition. Twenty-four such 
images were acquired over 360° with 15° resolution, (e) Two metrics for vortex displacement: (i) Width at 30% of the vortex 
height, corresponding approximately to the radius of maximum lateral force between tip and vortex, shown in green; (ii) 
Displacement between the peaks of the forward and backward scans, shown in blue, (f) Polar plot showing no signihcant 
anisotropy using either displacement metric. Polar plot data is shown for a single vortex, but is representative of similar plots 
for 5 different vortices. 


Hq = 300 nm. We note that our choice of cutoff height 
Zcut within the range 200-300 nm does not change the 
final result. 

Alternatively, we may obtain Meissner response curves 
from the height series of vortex images in Figs. 3(a-e). 
For each image that we fit, there is a constant offset as 
a fit parameter. We attribute this constant offset to the 
Meissner response of the superconductor to the tip. Plot¬ 
ting this constant offset with respect to the image tip 
height gives the equivalent of a Meissner response curve. 
Fitting the same model, Eq. D4, to the data points ob¬ 
tained this way gives a consistent value of Xab = 300 ± 60 
nm, as shown in Fig. 4(b). 

Our MFM measurements allow us to extract Xab 
several independent ways, all assuming a truncated cone 
tip shape. First, our fits of Eq. C4 to the vortex images 
exemplified in Eigs. 3 gave Xab ^ 500 nm. However, the 
finite lateral extent of the tip typically leads to an over¬ 
estimation of Xab using this kind of analysis.Second, 
the best fit of the {dFz/dz)max data in Eig. 3(f) to Eq. 2 
gives shows Xab = 380 ± 50 nm. This method is typically 
less sensitive to the details of tip shape.Einally, the 
Meissner analysis presented here in Appendix D depends 
only on data acquired far from vortices, where there is 
no lateral spatial dependence, and consequently very lit¬ 
tle sensitivity to the tip shape. We therefore find the 
Meissner method to be the most reliable measure of Xab 


Appendix E: In-Plane Anisotropy 

In-plane anisotropy in Ee-SCs can arise from 
anisotropy of the crystal structure,defects,band 
structure,or superconducting gap.^^ 

Structural Anisotropy 

Early phase diagrams for the ‘fill’ family of Ee-SCs, 
including La-1111 (Refs. 60-62), Ce-1111 (Ref. 63), Pr¬ 
ill 1 (Ref. 64), Nd-1111 (Ref. 65), and Sm-1111 (Refs. 
62 and 66) show no overlap between orthorhombic struc¬ 
ture and superconductivity. However, in most families of 
Ee-SCs, the orthorhombic phase overlaps with the super¬ 
conducting phase,and several more recent papers have 
questioned the non-overlap in Ce-1111 (Ref. 68) and Sm- 
1111 (Refs. 69 and 70), suggesting that weak orthorhom- 
bicity may persist nanoscale domains, making it hard 
to detect using standard procedures. The possible co¬ 
existence of orthorhombic phase and superconductivity 
in ‘1111’ materials remains controversial.Our own 
EBSD data in Eig. 1(d) cannot distinguish between or¬ 
thorhombic and tetragonal state in NdEeAsOi-a^E^,. 

Electronic Anisotropy 

Electronic anisotropy could arise from the underlying 
band structure or the superconducting gap. Theoreti¬ 
cal calculations predict negligible anisotropy of the bulk 
Eermi surface of NdEeAsOi-a^E^, (Ref. 73). This predic¬ 
tion is apparently born out in ARPES measurements,^^ 
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although one must be cautious because the ARPES- 
measured Fermi surface is significantly larger than would 
be expected for the bulk bands, and may not be repre¬ 
sentative of the bulk. ARPES also measured negligible 
superconducting gap anisotropy on NdFeAsOi_a,Fa,(Ref. 
33), although again one must caution that ARPES may 
be sensitive only to surface superconductivity which may 
not be representative of the bulk. There is some evidence 
that at least some Fe-SCs do have significant supercon¬ 
ducting gap anisotropy, such as the d-wave KFe 2 As 2 (Ref. 
75), but T-dependence of A in NdFeAsOi_a,Fa, seems to 
rule out the d-wave gap symmetry. 

Defect Anisotropy 

Fe site defects in Fe-SCs have shown anisotropy on two 
length scales. First, in both tetragonal and orthorhombic 
structural phases, Fe site defects form simple geometric 
dimers of two-unit-cell extent, randomly oriented along 
either of the two orthogonal Fe-As directions. These 
structures have been observed by STM in FeSe (Ref. 17), 
LiFeAs (Refs. 15 and 76), and LaFeAsO (Ref. 77). Sec¬ 


ond, in the spin density wave phase, Fe site defects form 
electronic dimers of size ^8 to ^IGape-Fe- These dimers 
are uniformly oriented along the a (longer, antiferromag¬ 
netic) axis. They have been observed directly in FeSe 
(Ref. 17) and NaFeAs (Ref. 78) and inferred indirectly 
in Ca(Fei_a,Coa,) 2 As 2 (Ref. 58). Their appearance is ex¬ 
plained by a competition between the impurity-induced 
(tt, tt) magnetism and the surrounding bulk spin density 
wave (tt, 0) phase. 

Search for Anisotropy in Vortex Wiggling 

To search for signatures of in-plane anisotropy, we in¬ 
vestigate the angle-dependence of vortex wiggling. This 
technique has been used to reveal pinning anisotropy in 
YBa 2 Cu 307 _a, (Ref. 51). Here we similarly raster the tip 
across a set of vortices at a series of angles, as illustrated 
in Figs. 9(a-d). The average wiggle distance of a vortex 
at each angle can be quantified in two different ways, as 
illustrated in Figs. 9(e-f). We found no statistically sig¬ 
nificant anisotropy of vortex wiggling using either metric, 
on 5 different vortices. 
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